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We employ ultra-broadband terahertz-midinfrared probe pulses to characterize the optical re¬ 
sponse of photoinduced charge-carrier plasmas in high-resistivity silicon in a reflection geometry, 
over a wide range of excitation densities (10^^ - 10^^ cm“^) at room temperature. In contrast to 
conventional terahertz spectroscopy studies, this enables one to directly cover the frequency range 
encompassing the resultant plasma frequencies. The intensity reflection spectra of the thermalized 
plasma, measured using sum-frequency (up-conversion) detection of the probe pulses, can be mod¬ 
eled well by a standard Drude model with a density-dependent momentum scattering time of ~200 fs 
at low densities, reaching ~20 fs for densities of ~10^^ cm“^, where the increase of the scattering rate 
saturates. This behavior can be reproduced well with theoretical results based on the generalized 
Drude approach for the electron-hole scattering rate, where the saturation occurs due to phase-space 
restrictions as the plasma becomes degenerate. We also study the initial sub-picosecond temporal 
development of the Drude response, and discuss the observed rise in the scattering time in terms 
of initial charge-carrier relaxation, as well as the optical response of the photoexcited sample as 
predicted by finite-difference time-domain simulations. 


I. INTRODUCTION 

The dynamics of charge carriers in silicon has been 
studied intensely over previous decades, due to its im¬ 
portance in (opto-)electronic applications. Spectroscopy 
provides an important tool for fundamental investiga¬ 
tions, as it allows one to elucidate the intra-band charge- 
carrier dynamics via the frequency-dependent conductiv¬ 
ity spectrum. For carriers in doped Si, this response has 
been characterized over a very broad range of n-/p-type 
dopant carrier concentration (10^^ - 10^^ cm“^), employ¬ 
ing measurements spanning the terahertz (THz) [T] and 
infrared [2] frequency ranges, where a Drude-type spec¬ 
tral behavior is generally observed, with modifications 
due to, e.g., inter-valence band transitions and energy- 
dependent carrier relaxation rates. A number of reports 
have also been devoted to the study of photoexcited 
electron-hole (e-h) plasmas, not only because of their rel¬ 
evance for applications (optoelectronics [3], laser micro¬ 
machining [4], extreme surface nonlinear optics 0) but 
also because they provide the possibility to study ultra¬ 
fast relaxation/scattering processes vs. excitation den¬ 
sity Wex and energy hz^ex- Such pump-probe studies typ¬ 
ically use inter-band excitation (UV-near-IR) and probe 
the spectral response in the THz 0 0, mid-infrared 
(MIR) [8] or optical [SHE] range, depending on (and 
hence the Drude plasma frequency z^p). Other studies 
have estimated the time scales directly from time-domain 
data using an ultra-short optical probe pulse m or pho¬ 
toelectron detection M- Among these reports, vari¬ 
ous approaches are used to estimate the scattering time 


Corresponding author: roskos@physik.uni-frankfurt.de 


Ts = F"^ from the data. In reviewing the range of reported 
data (presented in this paper), we found that there is 
some disparity in the values and their dependence on A'ex, 
and at high excitation densities (A^ex ^ 10^^ cm“^) there 
still appears to be contention about the order of magni¬ 
tude for Tg. For certain experiments (e.g. those involving 
the charge-carrier plasma as a moving front for Doppler 
up-conversion mm), the precise magnitude of Tq sen¬ 
sitively affects the degree of absorption loss for probe fre¬ 
quencies (especially about z^p), and hence a quantitative 
determination Tg vs. A^ex is necessary. While the role of 
the various scattering mechanisms has been investigated 
theoretically mm, a detailed comparison between the¬ 
ory and experiment is also lacking. 

In the present paper, we apply pump-probe spec¬ 
troscopy on Si(100) at T = 293 K in a reflection geom¬ 
etry, using both ultra-broadband MIR-probe and con¬ 
ventional THz-probe pulses, to extract estimates of Tg 
(following initial energy relaxation) over the range A^ex = 
3T0^^-2T0^^ cm“^, by applying a standard Drude model 
to fit the intensity spectra (and accounting for the longi¬ 
tudinal excitation density profile in the analysis). A satu¬ 
ration in the increase of the scattering rate is observed for 
A^ex'^lO^^ cm“^, with rg~20 fs. This behavior can be re¬ 
produced well by theory which includes the e-h scattering 
rate (whose variation dominates the density dependence 
in this range) via a generalized Drude approach [191 121 ] • 
This supports the assertion that the saturation is due to 
the onset of phase-space restrictions as the plasma be¬ 
comes degenerate (and the density of vacant final-states 
reduces). While the present time resolution is limited to 
a few hundred femtoseconds, we also investigate the sub¬ 
picosecond evolution of the plasma response. A compar¬ 
ison to model results from one-dimentional (ID) finite- 
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FIG. 1. (a) Spectrogram of ultra-broadband THz-MIR 

continuum probe pulse (black curve indicates temporal first- 
moment vs. frequency, (b) Corresponding spectral 

marginal M{y). (c) Estimate of temporal electric field pro¬ 
file E(t), obtained using M(z/) and (assuming a carrier- 
envelope phase of 920 = 0, see text). 

difference time domain (FDTD) simulations shows that 
while one should be wary of artifacts in the fitted values 
of Tg during the initial rise of the signal, the data in¬ 
deed indicate that the scattering rate relaxes to a steady 
state-value during the first 1 ps following excitation. 

II. EXPERIMENTAL DETAILS 

The experimental system is based on a 1-kHz 
Ti:Al 2 Os-amplifier laser (Clark-MXR, CPA-2101, Aq = 
775 nm, Tf^hm = 150 fs). For the ultra-broadband THz- 
MIR probe pulses, we employ the emission from a two- 
color air-plasma (as reported previously EaESD, whose 
pump beam (energy 380 /iJ) is obtained by spectral 
broadening of the Chirped Pulse Application (CPA) laser 
pulses in an Ar-filled (2.5 bars) hollow-core fiber and sub¬ 
sequent recompression with a set of negative-dispersion 
mirrors. The pump beam was focused to generate the 
plasma with a planoconvex lens (/ = 200 mm) via a 
150 /im-thick crystal (32°-cut) to provide the 

second-harmonic beam. The ultra-broadband emission 
had a typical pulse energy of 100 nJ (directly after the 
plasma, as inferred from measurements with a calibrated 
pyroelectric detector) and spectral coverage extending to 
-150 THz (5000 cm“^, 0.62 eV). This beam is collimated 
by an off-axis paraboloidal mirror (OAPM, effective focal 
length /eff = 152.4 mm) after passing a 450 /im-thick Si 
wafer (to discard the optical pump beams). The optical- 
pump THz-MIR-probe reflection measurements were per¬ 
formed by deflecting the probe beam with a 450 /im-thick 


Si beamsplitter to the sample, where it was focused by 
an OAPM (/eff = 101.6 mm) at normal incidence, re¬ 
flected, and transmitted through the Si beamsplitter to 
the detection stage. The CPA sample-pump beam (with 
pulse energy up to -50 /iJ) was brought to a diameter of 
0.6 mm (FWHM) on the sample, using a non-collinear ge¬ 
ometry. For detection of the THz-MIR probe pulses, we 
employed sum-frequency (SF) generation (up-conversion) 
[24] in a 500 /im-thick (llO)-cut ZnTe crystal with 150- 
fs optical detection pulses, whereby the THz-MIR beam 
is collinearly focused with an OAPM (/eff = 101.6 mm). 
Following the SF crystal, a notch filter is used to discard 
the input optical light and the SF spectrum is measured 
by coupling free-space into a miniature spectrometer with 
a cooled CCD sensor (Ocean Optics QE65). Integration 
times of 50-100 ms are sufficient to acquire a SF spec¬ 
trum for each value of detection delay t. A measured 
reference spectrogram of these pulses (where v is 

the frequency offset from z/q = c/Aq) is shown in FiglHa). 
Here the beam was sent directly to the detection ^nd 
hence represents that which is incident on the sample for 
the measurements below). In order to estimate the in¬ 
tensity spectrum /(z^), we also calculate the frequency 
marginal M{u) = /dt5'(z2, t) as shown in Fig. IHb). Note 
that the spectrogram (and hence M{u)) is affected both 
by the phase-matching response (which for such thick 
crystals varies as P(z/, z/o)'^l/A/c(z/, z/q) [H]): as well as 
a convolution with the spectral optical detection spec¬ 
trum (with width Az 2 = 3.8 THz here). Nevertheless, 
for the calculation of reflectivity spectra, M(z/) can be 
used with reasonable accuracy, as long as the spectral 
resolution is included in subsequent analysis. One can 
see that the measured signal extends from the cut-off 
of the notch filter (-10 THz) up to -150 THz. The 
modulation in the spectrum is due to the absorption 
lines from ambient H 2 O (-50 and 110 THz) and CO 2 
(-70 THz) in the beam path. While the temporal width 
of the spectrogram is dictated by the optical detection 
pulses, one can still extract an estimate of the group de¬ 
lay Tg(z 2 ) via the temporal first moment vs. frequency, 
i.e. t{u) = /dtt-5'(z/, t)//dtp(z/, t), which is also included 
in Fig. IHa) (black line), and indicates a pulse duration 
of -100 fs, which is due to primarily to the dispersion of 
the 450 /im-Si wafer. While only approximate, to give 
an impression of the time-domain pulse, in Fig. 00) we 
plot a reconstructed field profile using the approxima¬ 
tion F^(z/)->yM(z/)e“*-^(note that SFG detection 
can only provide the relative spectral phase, and we as¬ 
sume an arbitrary value of (/?o = 0 for the carrier-envelope 
phase here to illustrate the temporal chirp of the carrier 
wave). 

The purity of the high-resistivity Si samples used (float 
zone, (100), thickness of 525 /im, p > 3000 (1 • cm, Crys- 
tec GmbH) was characterized by two different spectro¬ 
scopic methods. Firstly, shallow group HI and V im¬ 
purities were analyzed by low-temperature near-IR pho- 
toluminescence at T = 4.2 K using a Bruker Vertex80 
FTIR spectrometer with a photoluminescence module 
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B -2.9-10“ P -7.9-10“ C <5-10^® 

A1 <5-10“ As <5-10“ O <5-10“ 

Sb <5-10“ 


TABLE L Determined elemental impurity concentrations (in 
cm“^) in Si samples from low-temperature photoluminescence 
(4 K) and mid-IR absorption (300 K) measurements. 


and liquid-helium immersion cryostat. Quantification 
of impurities was carried out according to the SEMI 
MF1389 standard [25] using a suitable set of calibra¬ 
tion samples. Excitation with two different wavelengths 
(Ai = 532 nm, absorption depth aex~l A 2 = 976 nm, 
<^ex'"70 jam) gave similar results, indicating no signifi¬ 
cant gradient of the shallow impurity concentration. Sec¬ 
ondly, substitutional carbon and interstitial oxygen con¬ 
tent were analyzed according to the SEMI MF1391 and 
SEMI MF1188 standards with room-temperature MIR 
absorption spectroscopy using a Bruker VertexTO FTIR 
spectrometer. The concentrations (upper bounds) of 
III/V dopants and C/0 are listed in Table which in¬ 
dicate that the sample purity is not corrupted e.g. by 
significant compensation doping, compared to the exci¬ 
tation densities used here. 


III. RESULTS 

A. Excitation density dependence 

Figure l^a) shows the measured intensity reflection 
spectra Ryy^r)jR q{i') (i.e., relative to the reflection Rq 
of the unpumped sample, which is a flat curve with Rq = 
0.30) at a pump-probe delay of r = 1 ps and a range 
of front-face {z = 0) excitation densities Aex = 0.35-2.1 • 
10^^ cm“^. 

Each reflection spectrum was fitted assuming a stan¬ 
dard Drude model for the plasma conductivity a{uj,z) = 
e^A^ex(^)/^eff(1 + where rueff = 0.15mo is 

taken for the effective reduced e-h mass, and Nq^(z) = 
Aexe“^®^^ accounts for the ^-dependent longitudinal exci¬ 
tation. It is well established for the reflectivity of such an 
excitation profile (and probe wavelengths X = cju > 
that one should not use simply the Fresnel field reflec¬ 
tion coefficient rp (i.e. based on the front-face density) 
for quantitative analysis, rather one must account for 
the distributed reflection of (and losses within) the ex¬ 
citation profile [26l [27|. The solution for the total field 
reflectivity at normal incidence can still be expressed in 
closed-form, which upon inspection of the formulas in 
Ref. [26l can be expressed as r = (l-n')/(l-tn'), where 
n' = nb + \/Ae^X/ 3 (^), nb = ^/^ is the background 
(unpumped) refractive index, Aero(co’) = -ia{uj,0)/€QUJ 
is the front-face photoinduced change in the complex 
permittivity, and = /^+i(^)//^(^) is the ratio of 

modified Bessel functions with /3 = 2iujn}j{caQ^)~^ and 
^ = 2iuj\/ Ae^o (cQ^ex)”^ • During the fitting, we convolved 
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FIG. 2. (a) Experimental intensity reflection spectra 

R{ia) = M(z/)/Mo(z^) for different values of optically excited 
plasma density Aex (Aex = 775 nm) and a pump-probe delay of 
T = 1 ps, and fitted model curves assuming a standard Drude 
dispersion model and accounting for the finite-depth expo¬ 
nential excitation density profile (with Oex = 125cm"Q. The 
inset shows the corresponding group delay At(z/) and model 
curve using parameters from fitting intensity reflection, (b) 
Corresponding parameters Ts (o), Up (□), from model fits in 
(a) vs. Aex. Also included are fitted values of Ts (A) when 
using only the SF spectra at t = 0, for later discussion in con¬ 
nection with the data in Fig. The inset depicts schematic 
of experimental reflection geometry. 


the raw model spectra with the spectral detection re¬ 
sponse (width Ah' given above). 

The fitted model curves are included in Fig. |^a) and 
are generally in close agreement with the experimen¬ 
tal data. The corresponding fit parameters (tq and 
= ujpl27T where cjp = cr(0,0)/ebreo) are shown in 
Fig. §b). We firstly note that the values of corre¬ 
spond closely to a dependence (solid curve), which 

supports the idea that Aex F (where F is the exci¬ 
tation fluence) and shows that higher-order excitation 
effects (e.g. two-photon absorption m) do not play a 
role in these measurements. The values of the Drude 
scattering time Tg exhibit a moderate trend of increas¬ 
ing scattering rate T = with Aex, ranging from 
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Ts = 30.7 fs (A^ex = 0.35 • 10^^ cm-^) to n = 19.7 fs 
(A^ex = 2.1 • 10^^ cm“^), as discussed below. We note 
that the fit residuals in Fig. [^a) typically deviated be¬ 
yond the random noise levels in the data, i.e. a small 
systematic disagreement exists for each dataset, preclud¬ 
ing a standard fit-parameter error analysis [28]. In or¬ 
der to obtain confidence intervals for th e fitted values 
of Ts, we analyzed the rms-misfit O = y/'ZiR - Rmod)‘^ 
in the (z^p,rs)-neighborhood of the fitted parameters and 
give here the bounds for Tq for the elliptical region where 
O = \/20rnin- Note that while we fitted here only the 
intensity curves = |r(z/)p, for the field reflectivity 

response r(z^) one also expects a spectral phase corre¬ 
sponding to a small but measurable group delay ATg(z^). 
This can also be estimated from the spectrograms via 
At{u)= t{u) -to(z/), as described above. In the inset of 
Fig. I^a) we plot the experimental Ai{u) data and the 
model curve obtained using the intensity-fit parameters 
for the case Agx = 0.89 • 10^^ cm“^, which are seen to be 
in good agreement and further demonstrates that the ex¬ 
perimental data correspond well to the Drude/reflection 
model used. 

In order to put these measured values of Tg in perspec¬ 
tive, in Fig.j^we plot the values along with various previ¬ 
ously reported estimates from the literature (as indicated 
in caption) vs. Agx- In order to provide additional data 
for lower Aex, we also performed a set of optical-pump- 
THz-probe reflection measurements using a conventional 
THz time-domain spectroscopy (TDS) system covering 
the range 0.1-3 THz. The obtained values of Tg (for 
Aex = 3.3 • 10^^ - 3.3 • 10^^ cm“^ at r = +5 ps) are also 
included in Fig. (where we applied the same model as 
above, except here fitting the complex reflectivity data 
r{u)). The literature data were obtained using various 
techniques, some of which do not directly probe precisely 
the Drude scattering, i.e., current damping, rate, includ¬ 
ing reports in which the values are extracted directly from 
analysis and those in which estimates are cited that are 
found to be consistent e.g. with the results of auxiliary 
simulations (and as such may only indicate an approxi¬ 
mate value of Tg). 

For low excitation densities, our results give a value 
of Tg = 200 ±75 fs, which is close to cited values for the 
e-ph scattering time F^.p^ at T = 300 K of 240-260 fs 
iniiiiiM]. A previous pump-probe THz-TDS study 
(■, Aex = 3 • 10^^ cm' -3 m estimated a value of Tg = 
330 - 500 fs (r/27r = 2-3 THz), although this report con¬ 
centrated on simulating the response during the initial 
dynamics (discussed further below) and accounting for 
the finite pulse duration of the THz probe pulse, as op¬ 
posed to a rigorous analysis of the conductivity spectra 
of a quasi-static plasma at sufficiently large delay r. A 
more recent study using THz-TDS [7] at T = 30 K with 
(Aex = 400 nm. A) fitted the complex conductivity spec¬ 
tra a{iy) for r = ±10 ps as a function of A^ex- Given the 
low temperature, the higher Tg values (~2 ps) at low den¬ 
sities are reasonable, as e-ph scattering is strongly sup¬ 
pressed [32]). However, the observed trend with A^ex dif- 



FIG. 3. Experimental and theoretical (Drude) scattering 
time vs. Aex- Values from present work (Aex = 775 nm) using 
THz-TDS (•, see text) and mid-IR continuum probe with 
SFG-detection (o, data from[^b)). Literature sources (incl. 
Aex, at room temperature unless otherwise stated): Ref. [6] 
(■, 760 nm); Ref. [7] (A, 400 nm, T = 30 K); Ref. M (□, 
800 nm); Ref. [14] (☆, 610 nm); Ref. [13] (★, 800 nm); Ref. |30] 
(♦, 530 nm); Ref. [T2] (X, 625 nm); Ref. [TO] (+, 620 nm). 
Gurves correspond to theoretical calculations of the total scat¬ 
tering rate (solid) from e-h (dashed) [21], acoustic- (dotted) 
and optical-phonon (dash-dot) scattering [T7], respectively 
(see text). 

fers significantly from our results, as discussed further be¬ 
low. Our data are consistent with the value of Tg = 100 fs 
employed in Ref. [29](n) used to simulate experiments in¬ 
volving the reflection from a counter-propagating plasma 
boundary in the range Aex = 1.5 • lO^^-lO^^ cm“^ (al¬ 
though no error margin for this value was given). 

In the range Aex = 10^^ - 10^^ cm“^, two reports pro¬ 
vide estimates of carrier scattering time which are rea¬ 
sonably consistent with our data, although they do not 
specifically probe the current relaxation rate. In Ref. jTT] 
(☆, Aex~10^^ cm“^), optical two-photon photoemission 
measurements were used to deduce an initial time scale 
for momentum relaxation of Tg = 40 fs, whereas in Ref. M 
(★, Aex = 5.5 • 10^^ cm“^) a time scale of Tg = 32 ± 5 fs 
was determined from the decay of the coherent four- 
wave-mixing signal. In both of these reports, this time 
scale accounts specifically for the initial (elastic) mo¬ 
mentum reorientation and redistribution of the carriers, 
while the energy relaxation time (due to e-ph scattering) 
was determined through additional measurements to be 
A-ph~250 fs. Hence these studies include the e-h and e- 
e scattering contributions directly, the latter of which is 
only expected to contribute to the Drude current damp¬ 
ing via the conduction-band anisotropy [21]. We note 
that another study [8] of MIR reflectivity spectra (simi¬ 
lar to the MIR measurements here only with a point-wise 
probing with tunable source) for Aex~10^^ - 10^^ cm“^ 
used a nominal value of Tg = 100 fs to provide a set of ad¬ 
ditional simulated curves. However, as no attempt was 
made to extract Tg from the experimental data, we omit 
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this value in Fig. 

The remaining data at higher excitation density A^ex > 
10^^ cm“^ all involved using a single probe wavelength 
and varying A^ex for determining the Drude response. 
The value of Tg = 100 fs in Ref. |30] (♦, A^ex~0 • 10^^ cm“^), 
which tends to deviate significantly from the overall trend 
of the data, was cited in connection with detailed simula¬ 
tions of the delay-dependent reflectivity of a 20-ps probe 
pulse (A = 2.8 /im), although no estimate of the error was 
provided and perhaps should be taken to indicate only 
the order of magnitude of Tg. The extreme excitation den¬ 
sities (and hence z/p values) of Ref. [12] (X, A^ex'^lO^^ cm 
and Ref. [10] (+, A^ex~10^^ cm“^) required the use of a 
visible-range probe, with estimates Tg^l.l fs and ~0.3 fs 
obtained by fitting the Drude-type reflection curves, re¬ 
spectively. While the fitted curves around v = those 
reports are in reasonable agreement with the experimen¬ 
tal data (although in both reports there are deviations 
on the high-A^ex side), it is not clear whether such an ap¬ 
proach can be applied robustly to extract Tg. In simula¬ 
tions, we found that the extracted value of Tg from reflec¬ 
tivity spectra is particularly sensitive to any distortion 
of the curves, as it is essentially the curvature of R{i') 
about z^p which allows its determination. Given that a 
Drude analysis of data with a single probe wavelength 
necessitates the variation of A^ex (and hence Tg), distor¬ 
tions of the measured R{Nq^) curves are to be expected 
(compared to spectra R(z^) at constant A^ex), and we ex¬ 
ercise caution in considering these values. We note that 
in both reports, it is the initial sub-picosecond carrier 
plasma response that is probed - i.e. in Ref. [12] R(A^ex) 
is measured with a pump-probe delay of r = 150 fs with 
100-fs pulses, while in Ref. uni the self-induced reflec¬ 
tion changes are measured with 90-fs pulses. However, 
in a preceding experimental report m, the authors of 
Ref. [12] and colleagues extracted Tg from their data vs. 
pump-probe delay (for a fixed A^ex~10^^ cm“^), asserting 
that Tg 0.5 fs for the whole measured delay range out 
to r = 0.8 ps. Note that for those measurements, a single 
probe wavelength was used with an analysis of p- and 
s-reflectivity at oblique incidence to extract the complex 
refractive index and hence Drude parameters, as opposed 
to analysis of R(A^ex)-curves. To summarize, while some 
concern has been raised about the plausibility of these 
values [21] , the assertion in the literature that scattering 
times Tg^l fs are obtained at these high densities (even 
for time scales > 100 fs) has not yet been refuted. 

B. Theoretical predictions 

We now turn to theoretical predictions for the observed 
dependence on A^ex- Qualitatively, the e-phonon scat¬ 
tering rate Te-ph is expected to vary only weakly with 
A^ex for densities below -10^^ cm“^ [32], whereas (bi- 
molecular) e-h scattering should obey Te-h ^ A^ex for 
low densities. While like-charge scattering does not con¬ 
tribute to current damping for isotropic bands, for the 


anisotropic X-valleys in Si e-e-scattering can also con¬ 
tribute, although in the low-frequency limit this is pre¬ 
dicted to have only a very minor influence on the mea¬ 
sured scattering time m- The situation for increasing 
A^ex was closely examined theoretically [m m El] fol¬ 
lowing the sub-femtosecond experimental value found in 
Ref. Iini Using different approaches, they all predict that 
as the plasma becomes degenerate (A^ex ^ ^ex^? where 
= (2mde^T/fz^)^/^/37r^, [17]) Te-h should reach a 
maximum value and then begin to decrease with A^ex, 
due to the reduction of available final states for scatter¬ 
ing (phase-space restrictions in the quantum limit [H]). 
For electrons in the conduction band at T = 300 K one 
has A^ex^ = 3.4-10^^ cm“^, using the value rride = 0.32 - me 
for the (valley-degeneracy-scaled) DOS mass. Although 
the effective temperature T can be expected to rise with 
excitation fluence due to the excess excitation photon en¬ 
ergy (for hz/ex > = 1.12 eV), at T = 300 K this amounts 

to only some 10 K for hr'ex = 1.6 eV. Hence one expects 
that the e-h scattering dominates the density dependence 
of Fex: and that a saturation of the F oc dependence 
should be observed for A^ex~10^^-10^^ cm“^, as is indeed 
indicated in our experimental results. 

We carried out calculations of Fe-h vs. A^ex using the 
theoretical treatment in Ref. HD i e. based on a general¬ 
ized Drude model, where the e-h scattering is calculated 
from integration over the appropriate expression of the 
complex polarizability functions (ae,h(q 5 ^) to yield the 
relevant coefficients of mutual friction, which are then 
used to evaluate the complex conductivity cr(z^). Here 
we assume isotropic bands for simplicity (with effective 
masses rude = 0.32 • tUq and ruh = 0.52 • ttIq) and calculate 
the low-frequency limit for Fe-h (which is a reasonable 
approximation for the spectral range used in our experi¬ 
ments) using the finite-temperature (T = 300 K) expres¬ 
sion for Im{(ae,h} and zero-temperature expression for 
Re{(ae,h}- The calculated rates Fe-h vs. Fgx are shown in 
Fig. (dashed curve), which indeed exhibit a saturation 
of the Fe-h ^ Aex dependence in range mentioned above, 
reaching a minimum with Tg 20 fs at = 8 • 10^^ cm“^ 
above which the e-h scattering rate begins to drop again. 

In order to include the contribution of e-ph scattering, 
we employed the theoretical expressions in Ref. [17] (Eqs. 
7 and 8 for the optical and acoustic phonon scattering, 
respectively). The calculated curves are also shown in 
Fig. (using the following values: LO optical phonon 
deformation potential Di = 6 • 10^ eV cm“^ [33]; op¬ 
tical phonon energy ^Doph = 64.4 meV [34]; acoustic 
deformation potential Ei = 8.2 eV [35]; sound velocity 
Vs = 9 ‘ 10^ ms“^). As expected, the acoustic phonon 
scattering rate Fe-aph dominates in the low-density limit, 
with a value Tg 250 fs reproducing experimental val¬ 
ues well. Based on this model, Fe-aph increases gradually 

(with a Fex -dependence due to an increase of the elec¬ 
tron velocity) reaching a maximum at 10^^ cm“^, 

upon which it decays with A^^^ due to the additional on¬ 
set of phase-space restrictions (as per Fe-h above). The 
contribution of optical phonon scattering is marginal, in- 
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creasing only at very high densities where electrons in 
the distribution increasingly reach energies 

The combined scattering rate is also shown in Fig. 
(solid curve), and is seen to reproduce the trend and 
magnitude of our experimental data well. To our knowl¬ 
edge, such a correspondence between experiment and the¬ 
ory for the Drude scattering time for carrier plasmas 
in undoped Si vs. A^ex has not been demonstrated un¬ 
til now. This correspondence has an important impact 
for judging the overall trend in the experimental data 
in Fig. On the basis of the experimental data alone, 
one might consider that the very short reported values 
of Tg for A^ex >10^^ could be consistent with a F oc Nq^- 
dependence extending to these high densities (as has been 
previously asserted HOD, with the apparent saturation 
of the values from our measurements for A^ex'^lO^^ cm“^ 
possibly being due to unforeseen artifacts in the experi¬ 
mental method. However, the consistency in the present 
results seems to rule this out, at least for charge-carrier 
plasmas after any carrier-carrier thermalization on a sub- 
100-fs time scale (see next section). Clearly the treat¬ 
ment here of e-h and e-ph scattering does not predict 
scattering times of Tg^l fs at high density. We also have 
difficulties to reconcile the data from Ref. [7| (Fig. A). 
In considering the low temperature for those measure¬ 
ments, one expects the onset of plasma degeneracy at 
even lower densities than for our experiments. Further 
experiments vs. T and excitation energy are required to 
resolve this issue, where a broader probe spectral range 
including the plasma frequency is covered, as per our 
measurements here. 


C. Sub-picosecond relaxation dynamics 

In this last section, we present time-resolved measure¬ 
ments of the plasma response probed with the THz-MIR 
pulses. As will be shown, while certain aspects of the 
data would necessitate superior time-resolution than that 
available here, we can draw important conclusions re¬ 
garding the analysis and interpretation of such measure¬ 
ments via comparison with simulations. 

The time-resolved THz-MIR reflectivity spectra 
R'{iy,r) were measured in the same experimental geome¬ 
try as the data in Fig. I^a) with A^ex = 2.1-10^^ cm only 
here, in order to achieve practical measurement times 
we acquired each SF spectrum at a fixed time 

t = to, as opposed to acquiring a full spectrogram to cal¬ 
culate the marginal (hence the measurements capture the 
spectral components within the 150-fs SF detection gate 
pulse). The data are shown in Fig. ^a). For each spec¬ 
trum, we employed the same Drude model analysis as 
above, in order to yield and Tg as a function of r. 
The model spectra and corresponding fit parameters are 
shown in Fig. [^b,c), respectively. As can be seen, the 
data can be reproduced well over the full delay range with 
this model. The kinetics of the fitted A^ex-data can be fit 
well assuming a Gaussian response (i.e. error-function 


dependence) with a FWHM width of Tr = 420 fs (due to 
the combination of the pump-probe correlation function 
and non-collinear geometry). The corresponding kinetics 
of the Tg-values show a rise from rg~17 fs around r = 0 to 
reach a value of Tg = 28 fs for a delay of some ps. The 
continuing rise of the Tg-curve compared to that of A^ex 
clearly indicates a subsequent evolution of the plasma 
scattering rate following excitation. This trend is not 
consistent with that expected simply for an increasing 
plasma density during the excitation pulse, and is ad¬ 
dressed in detail below. A comparison of the value for 
r = +1 ps (rg 26 fs) with the corresponding data in 
Fig. ^b) (Tg 19.7 fs) reveals an apparent discrepancy 
(both sets of data were acquired in the same measure¬ 
ment run). As the only difference in the two measure¬ 
ments is that in one case the full SF marginal M{u) is 
acquired vs. the SF spectra at fixed t = to, as a first step 
to address this issue we re-analyzed the data in Fig. 
by extracting only the single spectra 5'(z/, to) from the 
full spectrogram. The resultant fitted values of Tg are 
also included in Fig. |^b) (solid triangles), and are in¬ 
deed systematically larger than those obtained with the 
full spectrogram marginals (and then are consistent with 
the value in Fig. [^c) for Atgx = 2.1-10^^ cm ^). Hence an 
important finding is that SF time-gating introduces ar¬ 
tifacts in such measurements; we explain the reason for 
this below with the help of simulation results. 

We modeled the experiments using ID FDTD simula¬ 
tions, with a Drude response and spatio-temporal pro¬ 
file for the photoinduced plasma (as described previously 
na na, except here with co-propagating pump-probe 
pulses and front-face excitation). Simulation batches 
were carried out vs. the pump-probe delay r, using a 
time-independent nominal scattering time of Tgo = 25 fs. 
The numerical reflected fields were then analyzed to pro¬ 
duce either the full spectrum or that corresponding to 
SF detection at a fixed time to (with a 150-fs SF de¬ 
tection gate), corresponding to the various experimental 
results here. These spectra were then fitted using the 
same Drude model algorithm. We present the results for 
three sets of simulation conditions in order to isolate var¬ 
ious aspects expected in the analysis of the experimental 
data, as shown in Fig. [^d) (curves (i-ih)). 

In the first case (i), we used a transform-limited (TL) 
THz-MIR probe pulse with an intensity spectrum based 
on the experimental one (Fig. [^, a pump pulse dura¬ 
tion equal to the experimentally observed for A^ex in 
Fig. I^c), and fitted the full spectra of the numerical re¬ 
flected THz-MIR pulses for each r. The kinetics of the 
fitted A^ex-values are as expected and correspond closely 
to the experimental data. Around zero-delay, the fitted 
values of Tg are also seen to rise from < 15 fs toward 
the simulation value of Tgo = 25 fs (note the slight error 
is due to the finite accuracy of the FDTD simulations). 
However, it is also apparent that this artificial rise fol¬ 
lows the same delay-dependence as the Aex-curve, i.e. 
the effect occurs only during the formation of the plasma 
when there is temporal overlap between pump and probe 
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FIG. 4. (a) Experimental relative intensity reflection spectra 

vs. pump-probe delay r, estimated using SF spectra at fixed 
detection time to, i-e. to;7 ")/aS'o(z^, to), where 

is the reference spectrum before excitation (measured 
at sufficiently negative r). (b) Corresponding model spec¬ 

tra using Drude model fits for each value of r. (c) Fitted 
Drude parameters Ts (o), iVex (□) for model data in (b). (d) 
Corresponding analysis of FDTD simulation results with a 
constant value of Tso = 25 fs, demonstrating artifact rise of fit¬ 
ted Ts values around zero delay for three simulation batches: 
(i) transform-limited (TL) THz-MIR pulse with full-spectra 
(solid), (ii) including experimental THz-MIR pulse chirp and 
simulating SF detection at fixed delay (dashed), and (hi) with 
a TL probe pulse and a sequence of static, partially formed 
plasma profiles (short dash). 


pulses. 

The second case (ii) adds the additional effects of our 
measurements in Fig. ®a-c), by adding the experimental 
chirp = 126 fs^) to the input THz-MIR pulse, and 
using the numerically time-gated SF spectra. This leads 
again to a rise in the fitted Tg-curves, albeit with a more 
complex variation during the leading portion of the ex¬ 
citation pulse, r < 0. (Note the fitted A^ex-curve is essen¬ 


tially the same for all simulation cases, and only that for 
case (i) is shown). However, in this case, one sees that fit¬ 
ted Tg values in the steady-state are systematically larger 
than the underlying value of Tgo used in the simulations 
(as seen in the experimental analysis above). This is not 
the case when the full spectra are used (with either TL or 
chirped THz-MIR pulses, latter results not shown). This 
allows us to identify the origin of the discrepancy above: 
the effect of the finite SF-detection gate pulse introduces 
a small distortion on the measured spectra which results 
in a systematic error (overestimation) in the fitted value 
of Tg, despite that fact that the gate is also applied to 
the detection of reference pulse (without plasma excita¬ 
tion). Having the detailed simulation results at hand, 
this can be traced to the small group-delay variation ex¬ 
perienced by the THz-MIR pulse in reflection (see inset 
in Fig. [^a)). While this is only of the order of ±10 fs, 
it is indeed sufficient to distort the spectra in this way. 
This is an important result, as one may have expected 
that the use of a reflection geometry should remove any 
influence of plasma dispersion (compared to bulk delays 
in transmission), and emphasizes the role of the finite 
group-delay in reflection from (and associated penetra¬ 
tion into) a Drude medium. In future investigations this 
effect could be alleviated by the use of a significantly 
longer SF detection pulse (i.e. Tsf ^ 1 ps via spectral 
filtering). Nevertheless, one sees that rise of the fitted 
Tg-curve for case (ii) still is essentially complete after the 
pump pulse excitation, and that these effects alone do 
not account for the slower rise in the experimental data 

(Fig. Be)). 

Now that we have demonstrated numerically that such 
experiments inherently produce artifacts in the fitted 
value of Tg around zero-delay (which could cause phys¬ 
ical misinterpretation of data, even with superior time- 
resolution), the question of its specific cause remains. 
About zero-delay there are two deviations from the 
steady-state situation assumed in the Drude model here 
(besides the possibility of a time-dependent Tg-value), 
i.e. not only does the THz-MIR pulse co-propagate with 
the pump pulse and experience a time-nonstationary 
medium, but also the excitation profile assumed in the 
model (exponential decay with z) is not yet fully formed 
and hence the expression used for calculating r{i') (see 
Sec. lIlL^ [26]) is not completely accurate. In order to 
isolate these two factors, we also performed a simulation 
batch (Fig.j^d), case (iii)) with a control situation where 
the plasma for each value of delay is kept artificially 
“frozen” at its profile for t = r (and the THz-MIR pulse 
reflects from a t\m.e-stationary^ partially formed plasma 
profile). As can be seen from the fitted Tg-curve (iii), 
while one still observes an artifact rise (due to the inac¬ 
curacy of the plasma profile model), this is now signifi¬ 
cantly more abrupt. Hence one may conclude that the co¬ 
propagation and time-nonstationary medium dominate 
the slower rise seen in the curve for case (i). This is an 
important result for interpretation of future experiments, 
as this effect will be present even for much shorter pump 















































and probe pulses, and one should also employ detailed 
simulations of the co-propagation to avoid misinterpre¬ 
tation of experimental data. 

Finally, returning to the delay-dependence of Tg in 
Fig- 0c), in Fig. B c) we provide a fit of the data (for r > 
200 fs) with an exponential step function oc (1 - 
which yields a value of Tg = 450 fs. Given the simulation 
results above, it is clearly difficult to define and decon¬ 
volve a response function from the Tg-curve, and attempts 
to fit the Tg-curve including convolution with a Gaussian 
response degraded the quality of the fit. In comparing 
this value of Tg with the fitted Gaussian response for 
A^ex 5 one should consider the half-width-half-maximum 
of Tr/2 = 210 fs. Still, this uncorrected value of Tg is 
close to Tg = 500 fs reported in Ref. [6| (i.e. for a lower 
density, A^ex = 3-10^^ cm“^), and quite likely corresponds 
to subsequent energy relaxation of the photoexcited car¬ 
riers. 


IV. CONCLUSIONS 

The use of THz-MIR spectroscopy with a detailed 
quantitative analysis has yielded a more comprehensive 
estimate of the Drude scattering time of charge-carrier 


plasmas in undoped Si and its dependence on The 
momentum relaxation time of the thermalized plasma 
ranges from ~200 fs at low electron-hole densities to 
-20 fs for Vex~10^^ cm“^, at room temperature. The 
direct comparison with theoretical predictions provides 
mutual support for both, and supports the mechanism 
of scattering rate saturation due to phase-space restric¬ 
tions at high density. A review of the literature reports 
indicates that there are still open issues to address, both 
experimentally and theoretically. In particular, (i) the 
effect of excitation energy and temperature (in order 
to account for the differences between our results and 
those from e.g. Ref. o, (ii) the magnitude of Tg at very 
high densities (where either a revision of the experimen¬ 
tal results, or additional physical scattering mechanisms, 
would be needed to reconcile the reported values Tg-l fs), 
and (hi) the precise sub-picosecond plasma dynamics. In 
the latter case, state-of-the-art theoretical methods to 
describe the nonequilibrium dynamics [36] are required. 
Our ongoing THz-MIR studies with improved time reso¬ 
lution are aimed at providing an experimental test for 
such theoretical predictions. The results from auxil¬ 
iary ID-FDTD simulations revealed important aspects 
for conducting a reliable analysis of time-resolved spec¬ 
tra, and demonstrate that the spatio-temporal propaga¬ 
tion effects must be taken into account even in a reflection 
geometry. 
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